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Abstract. 

In this work we study the microwave photoconductivity of a two-dimensional 
electron system (2DES) in the presence of a magnetic field and a two-dimensional 
modulation (2D). The model includes the microwave and Landau contributions in 

£N| , a non-perturbative exact way, the periodic potential is treated perturbatively. The 

Landau-Floquet states provide a convenient base with respect to which the lattice 
potential becomes time-dependent, inducing transitions between the Landau-Floquet 

f^ | levels. Based on this formalism, we provide a Kubo-like formula that takes into 

account the oscillatory Floquet structure of the problem. The total longitudinal 
conductivity and resistivity exhibit strong oscillations, determined by e = uj/lo c with 
to the radiation frequency and lo c the cyclotron frequency. The oscillations follow a 

i-rt ' pattern with minima centered at u>/uj c = j + i(l — 1) + 5, and maxima centered at 

C | lo/lo c = j + \{l — 1) — 5, where j = 1,2,3 , 5 ~ 1/5 is a constant shift and I 

is the dominant multipole contribution. Negative resistance states (NRS) develop as 
the electron mobility and the intensity of the microwave power are increased. These 
NRS appear in a narrow window region of values of the lattice parameter (a), around 

j^ ■ a ~ Ib, where lg is the magnetic length. It is proposed that these phenomena may 

be observed in artificially fabricated arrays of periodic scatterers at the interface of 
ultraclean GaAs / Al x Ga\- x As heterostructures. 



PACS numbers: 72.40,+w, 73.21.Cd, 75.47.-m, 73.43.-f, 

1. Introduction. 

The microwave irradiation of two-dimensional electron systems (2DES) has remarkable 
consequences on the transport properties at low magnetic fields. Recently, two 
experimental groups [H El El El , reported the unexpected discovery of zero resistance 
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states (ZRS) when high mobility GaAs / Al x Gai_ x As heterostructures in weak magnetic 
fields were exposed to millimeter irradiation. Unlike the strong magnetic field regime, 
the Hall resistance is not quantized. The magnetoresistance exhibits giant oscillations, 
periodic in e = uo/uo c with uj the radiation frequency and uo c the cyclotron frequency; 
the series of minima formed at e = j + 5, j = 1,2,3...., 5 = § [XJ E], or £ = | 
H] . These discoveries triggered a large number of experimental 13 El E] and 
theoretical [TDl HH H21 [13 H31 [T31 [THl [TT1 UHl CHI I2DI I2U studies. According to reference 
[T5] . ZRS probably originate from negative resistance states (NRS); it was argued that 
negative resistance induces the formation of current domains, yielding an instability 
that drives the system into a ZRS. The existence of NRS was first predicted in the 
pioneering work of Ryzhii [TU1 ITT] . Nowadays two distinct mechanisms that produce 
negative longitudinal conductance are known: (i) the impurity scattering mechanism, 
which is caused by the disorder assisted absorption and emission of microwaves 
[JUJITTJ H2 EH H3J UJ)J Hnj[T71[THJ, and (ii) the distribution function mechanism, according 
to which the microwave absorption modifies the electron distribution function leading 
to a negative longitudinal conductance HHJ |20J l2T] . A model for the impurity 
scattering mechanism was proposed previously by the authors Jd QS], the model is 
based on the fact that the microwave and Landau dynamics can be exactly taken into 
account producing well defined Floquet-Landau states. The disorder effects are treated 
perturbatively, inducing transitions between the Floquet-Landau levels. The model 
reproduce various of the experimentally observed features, in particular the fact that 
negative resistance states (ZRS) appear only when the electron mobility exceeds a 
threshold. 

Although the experiments described above do not include the effect of periodical 
potential modulations, exploring its physical consequences is worthwhile. We can 
identified at least three reasons for doing so: (i) The study of both weak [22] and 
strong J2H] periodically modulated 2DES in the presence of magnetic fields has lead to 
the discovery of interesting transport effects, such as commensurability phenomena and 
transport anisotropies. (ii) The theoretical method previously developed in references 
[T7| ITH] is well suited to study this kind of system. (Hi) The use of artificially 
fabricated arrays of periodic scatterers at the interface of ultraclean GaAs / Al x Ga\- x As 
heterostructures |22l EH 123 EH] may allow to test the predictions made by these 
theoretical studies. 

In this work we make a theoretical study of the microwave photoconductivity 
of a 2DES in the presence of a magnetic field and a two-dimensional modulation. 
Two theoretical studies of modulated 2DES under the combined effects of magnetic 
and microwave radiation have recently appeared. Dietel et.al. J2H] considered the 
photoconductivity in the case of ID periodic modulation. The calculation uses first 
order perturbation theory for both the microwave field as well as for the periodic 
potential. Additionally, the calculation is restricted by the following conditions: the 
lattice parameter a is small as compared to the cyclotron radius R c , the temperature T 
is large with respect to the periodic strength potential V , and V <C tko c . Due of the 
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unidirectional structure of the modulation the photocurrents parallel and perpendicular 
to the modulation are different. The work of Gumbs [27] applies for a strong 2D 
modulation, but is linear with respect to the microwave field intensity. His approach 
made use of the usual Kubo formula in which the matrix elements are evaluated using 
the numerically obtained Hofstadter-type wave functions. In this paper we address 
the case of a 2D modulation under different conditions; in particular we explore the 
possible appearance of negative resistance states. The Landau and microwave field 
contributions are exactly taken into account. We analyze the regime in which the 
following conditions hold: kT ~ Vo < fajj C) anduT tr ~ uj c r tr » 1. Instead of appealing 
to the usual Kubo formula, our approach shows that the use of the Floquet-Landau 
states determines not only the wave functions required to evaluate the matrix elements; 
but also leads to a modified Kubo-like formula, in which the oscillations on the density 
of states induce the development of negative resistance states. As a first step we find a 
unitary transformation that exactly takes into account the dynamics associated with the 
Landau and radiation contributions. As a second step, the periodic potential is added 
perturbatively. With respect to the Landau-Floquet states, the periodic potential act as 
a coherent oscillating field which induces transitions between these levels. Based on this 
formalism, we provide a Kubo-like expression for the conductance that incorporates the 
oscillatory Floquet structure of the system. It is found that both a xx and p xx exhibit 
strong oscillations determined by e = uj/uj c . NRS develop for sufficiently high electron 
mobility and strong microwave power. The model is used to test chirality effects induced 
by the magnetic field, calculations are carried out for various E-field polarization's. 
Finally, we explore the nonlinear regime in which multiple photon exchange play an 
essential role, as well as the current-voltage characteristics of the system. 

The paper is organized as follows. In the next section we present the model and the 
method that allow us to obtain the exact solution of the Landau-microwave system, as 
well as the perturbative corrections induced by the periodic potential. In section El we 
develop the formulation of the linear response theory valid in arbitrary magnetic and 
microwave fields. A discussion of relevant numerical calculations is presented in section 
(J3J). The last section contains a summary of our main results. 

2. The Model. 

Let us consider the motion of an electron in two dimensions subject to a uniform 
magnetic field B perpendicular to the plane and a constant electric field E c , a periodic 
potential V and driven by microwave radiation. On the plane the dynamics is governed 
by the Schrodinger equation 

ih— = H^= [H {BM + V(r)] * . (1) 

Here H{b,u)} is written in term of the covariant derivative 

h { bm = ^n 2 , n = p + eA , (2) 
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where m* is the effective electron mass over the plane that takes into account the effects 
of the crystalline atomic structure over the charge carriers. The vector potential A 
includes all the contributions of the magnetic, electric and radiation fields: 



A = --r xB + Re 

2 



w exp{— iuii} 



E c £. (3) 



u> 
The superlattice potential V(r) is decomposed in a Fourier expansion 

V(r) = J2 V ™ exp {i27r (^ + M) } . (4) 

ran 

We first consider the exact solution of the microwave driven Landau problem, the 
periodic potential effects are lately added perturbatively. Along this work we shall 
assume: {%) a weak modulation | V^| <^C hu c and (ii) the clean limit UT tr ~ u c T tr » 1; 
here r ir is the transport relaxation time that is estimated using its relation to the electron 
mobility // = eT tr /m*. Based on these conditions the use of first order perturbation 
theory in V becomes reasonable. Furthermore the second conditions justifies the use 
of a Fermi distribution function evaluated at the Floquet-Landau quasi-energies, see 
appendix (J5.2J) . 

The system posed by H{b,lj} can be recast as a forced harmonic oscillator, a problem 
that was solved long time ago by Husimi |2HJ - Following the formalism developed in 
references "2~"H I3T"| . we introduce a canonical transformation to new variables Q u ,Pa] 
/i = 0, 1, 2, according to 

Qo —t, Pq = id t + e<f) + er ■ E, 

VeBQ 1 = U y , \^BP 1 = U X , 

Ve~BQ 2 = U X + eBy , VeBP 2 = U- y - eBx. (5) 

It is easily verified that the transformation is indeed canonical, the new variables obey 
the commutation rules: —[Qo,Po] = [Qi,P\] = [^2,^2] = iB; all other commutators 
being zero. The inverse transformation gives x = Ib (Qi — P2), and y = Ib (Q2 ~ Pi), 
where Ib = \ -^ is the magnetic length. The operators (Q2, P2) can be identified with 
the generators of the electric-magnetic translation symmetries pH] 13*2"] . Final results 
are independent of the selected gauge. From the operators in Eq. (J""} we construct two 
pairs of harmonic oscillator-like ladder operators: (ai,aj), and (02,02) w hh: 

01 = -i= (Pi - iQi) , a 2 = -j= (P 2 - iQ-2) , (6) 

obeying: [01, a\] = [02,02] = 1, and [01,02] = [oi,a 2 ] = 0. 

It is now possible to find a unitary transformation that exactly diagonalizes H^b,w}, 
it yields 

W^H {BM W = u c Q + a\ a^j = H , (7) 

with the cyclotron frequency u c = eB/m* and the W(t) operator given by 

W(i) = exp{ir]iQi} expli^Px} exp{ir) 2 Q 2 } exp{i£ 2 P2} exp{i / Cdt'} , (8) 
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where the functions rji(t) and £,(£) represent the solutions to the classical equations of 
motion that follow from the variation of the Lagrangian 



Ulr 



£ = y M + Ci) + Cim + C 2 ^2 + d B [E x (Ci + m) + E y ( Vl + c 2 )] • (9) 

It is straightforward to obtain the solutions to the equation of motion, using the 
expression for the electric field E = —dA/dt with A given in (J3J). Adding a damping 
term that takes into account the radiative decay of the quasiparticle, they read 



771 = el B E^ Re 
Ci = cIbEu Re 



-iue x + u c e y &iut 
ui 2 -ui 2 + iuT rad 
uj c e x + iuje y ' 



ui 



UI 2 + iuTrad 



rj 2 = el B E w Re 
C2 = -el B E w Re 



e y e 



iuit 



lul 



+ d B EH, 



e x e 



iujt 



lul 



-el B E c x t. (10) 



According to the Floquet theorem, the wave function can be written as *(£) = 
exp (—iSfjt) 4>n(t), where M (t) is periodic in time, i.e. ^(t+r^,) = M (t), with t w = 2n /ui. 
From Eq. (jHj) it is noticed that the transformed wave function ty w = W^f contains the 
phase factor exp f i J Cdt' J . It then follows that the quasienergies and the Floquet 
modes can be deduced if we add and subtract to this exponential a term of the form 
— /J" Cdt'. Hence, the quasienergies can be readily read off 

_ e 2 E 2 Jl + 2uj c Re(e x e y )/uj] 



£ u = £<?>+£ 



rad 



4°) = ^ c (-+/i 



2m* [( W - uj c ) 2 + T 2 rad ] 

here £l are the usual Landau energies, and the induced Floquet energy shift is given 
by the microwave energy £ ra d- The corresponding time-periodic Floquet modes in the 
(Pi,P 2 ) representation are given by 

* M , fe (P) = exp{-* sin (2cot) F(uj)} ( J) im (P 1 )S(P 2 - k) , (12) 

the index k labels the degeneracy of the Landau-Floquet states, and M (Pi) is the 
harmonic oscillator function in the P\ representation 

1 



{Pi) = (PM 



l l 2 2^\ 



e-^ 2 ^(Pi) 



T\ 



# M (Pi) is the Hermite polynomial and the function F(ui) is given as 



Ulr 



Fui) = - 



eEJ 



B 



id \ UJ 2 — LU 2 



UJ 



ui+2ujx-^y y + 



Re(e x e y 



ijJOJr 



(2lu 4 



2 2 
UI Ul„ 



+ - C 4 ) 



(13) 



.(14) 



Let us now consider the complete Hamiltonian including the contribution from 
the periodic potential. When the transformation induced by W(t) is applied, the 
Schrodinger equation in Jl} becomes 

p *W = #0*0*0 + V w (t)^ w) , (15) 

where * (w) = W(t)V and V w {t) = W(t)V(r)W~ l (t). Notice that the periodic potential 
acquires a time dependence brought by the W(t) transformation. The problem is now 
solved in the interaction representation using first order time dependent perturbation 
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theory. In the interaction representation W t = exp-fif/oO 1 ^ > an< ^ the Schrodinger 
equation becomes 

id t ¥^^{V w (t)}^fK (16) 

The equation has the solution Wj '(t) = U(t — t )Wj (to), where U(t) is the evolution 
operator. To first order in perturbation theory it is given by the expression 

U(t) = 1 - i f dt! [W\f)V(v)W(t , )\ I . (17) 

J — oo 

The interaction is adiabatically turned off as to ~~ * — oo, in which case the asymptotic 
states are selected as the Landau-Floquet eigenvalues of H , i.e. \^f\ (to)) — > \n,k). 
Utilizing the explicit expression for the W transformation in (JHJ) and after a lengthly 
calculation the matrix elements of the evolution operator can be worked out as 

f il B q^ ) (k+l B q' n v) /2) i(£^+ujl)t 

(/*, k\ U(t) \v, k') = s,j kk ,-J2 J2H k ~ k ' + rf) P - r- cj2 >fnB , (is) 

I rnn ^ 



where uje = el%(Qn E° — qm Ey), and the explicit expression for C^l^mn is given by 



CgU. = l 2 B V mn D, v (q mn ) ( ^^ ) J L (|A mn |) . (19) 



x mn\ 



In the previous expressions the discreet pseudomomentum are given as 

q^=2mn/a, q^ = 2nn/b, 

qmn^ih^-iq^/Vz, (20) 

Ji denote the Legendre polynomials and D^ v is given in terms of the generalized Laguerre 
polynomials according to 



1 U|2 



Q*T~\ %Lt-» (\q\ 2 ) , ul> 



v- 



D»»(q) = (v\D(q)\n) = e-m { V » " ' (21) 

and 

ui(l B eEu 



, » , - r ^K^z + g^eJ+Mg^-g^eJ] . (22) 

Summarizing, the solution to the original Schrodinger equation in Eq. (JI} has been 
achieved by means os three successive transformations: 

|*mW) = ^ exp{-iH t} U(t - t ) |/i, k), (23) 

the explicitly expressions for H , W, and £/ are given in Eqs. (JTj), (JHJ), and (fTSj) 
respectively. 
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3. Kubo formula for Floquet states. 

The usual Kubo formula for the conductivity must be modified in order to include 
the Floquet dynamics. In the presence of an additional DC electric field the complete 
Hamiltonian is Ht = H + V ex t, where H is the Hamiltonian in Eq. (£Q) and V ex t = 
—Tl-A ext , with A ext = — sin (fit) exp (—i]\t\) . The static limit is obtained with fi — > 0, 
and n represents the rate at which the perturbation is turned on and off. In order to 
calculate the expectation value of the current density, we need the density matrix p{t) 
which obeys the von Neumann equation 

ih^ = [H T ,p] = [H + V ext ,p]. (24) 

We write to first order p = p + Ap, where the leading term satisfies the equation 

iH^ = [H 7Po \. (25) 

In agreement with Eq. (J2BJ), Ap is transformed to 

Ap(t) = U\(t - t ) exp{iH t}W(t) Ap(t) W\t) exp{-zif t}^/(t - *o) • (26) 

In terms of the transformed density matrix Ap(£), Eq. (J24"jl becomes 
_dAp 

where V ext and p~o are the external potential and quasi-equilibrium density matrix 
transformed in the same manner as Ap in Eq. (|26|) . The transformed quasi-equilibrium 
density matrix is assumed to have the form p = ^ |p)/(^)(p|, where f{£^) is the 
usual Fermi function and S^ the Landau-Floquet levels. The justification for selecting 
a Fermi-Dirac distribution in the quasi-energy states is presented in the appendix (|5.2j) . 
If we consider the region in which the conditions r^ <C r tr <C r in hold, then the elastic 
and inelastic relaxation processes can be neglected as compared to the microwave field 
effects. The solution of the Boltzmann equation yields, a Fermi-Dirac distribution in 
the quasi-energy states [T7j, see appendix (|5.2|) . It is straightforward to verify that this 
selection guarantee that the quasi-equilibrium condition in (|25jl is verified. Using the 
results in Eqs. (|2*3*|l and (|2Tj). the expectation value of the density matrix can now be 
easily obtained from the integration of Eq. (|2"?j) with the initial condition Ap(£) — ► as 
t — > — oo giving for t < 

(tt M | Ap(i) 1*^) = (p, k\ Ap(t) 1 1/, k') 

eE, 



ih- 



V ext ,Po , (27) 



2 



„i(Cl—irf)t > 

f^^ k \U(t')\^ u ^) + (Q ^ -Q) 



(2f 



where the definition /^ = /(f^) — /(£,/) was used. The expectation value for the 
momentum operator is explicitly computed with the help of Eqs (jHJ), (fTSj) . and (|23*|) . 
retaining terms linear in the modulation potential and after a lengthly calculations it 
yields 

(* Mfe | II, I*,*,) = v 7 ^ J] J2 5 ( k ~ k ' + 1b ^ ] ) e"*^™^™ e ^ + " 1 -^ A% mn (j).(29) 
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Here the following definitions were introduced: £^ v — £ M — £ u , a, = bj = 1 if j = x, 

a i = ~ b j = ~i if J = V' and ^l,mn(j) is given by 



A« (7) = -— 



Sav — lo c + col — in £„ v + u c + u>l — in 



(30) 



the expression for C\w >mn is given in (JT§j) . It should be pointed out that in principle 
there is a zero order contribution (independent of V) to the expectation value of the 
momentum operator in ()29j) ; this would contribute to the direct cyclotron resonance 
heating that has a single peak around uj ~ u> c . However this contribution has proved to 
be negligible ^7], so it will be altogether omitted. Utilizing these results the time 
integral in Eq. ()28j) is readily carried out. The current density to first order in 
the external electric field can now be calculated from (J(£,r)) = Tr Ap(t)J(t) , the 
resulting expression represents the local density current. Here we are concerned with 
the macroscopic conductivity tensor that relates the spatially and time averaged current 
density j = (t w A)~ f " dt J d 2 x(J(t, r)) to the averaged electric field; here A is the area 
of the system (it is understood that A — > 00). Assuming that the external electric 
field points along the x-axis the macroscopic conductivity can be worked out. The total 
conductivity is given by a sum o x { = a® + o~ xi ; the dark conductivity is calculated 
when both the modulation and microwave radiation are switched-off as 



o 



i> e\ 2 r f/F 



xi ' Ahi ^ I fi 



in' 



(HfiS^u+i M<W-i 



Eya, + Q - in £y V + Q - in_ 



(il - -flU , (31) 



whereas the microwave-modulation (MM) induced conductivity is worked out as 
(MM) _ g u c V^ J ffW \^ V^ A/^ran (i)A yftmn (x) I 

)jlv \ ran I h ) 

Selecting i = x or i = y the longitudinal and Hall conductivities can be selected. The 
denominators on the R.H.S. of the previous equations can be related to the advanced 
and retarded Green's functions G^(£) = 1/ (£ — £^ ± in). To make further progress the 
real and absorptive parts of the Green's functions are separated taking the limit n — > 
and using lim v ^ l/{£ — in) = Pl/£ + iir5(£), where P indicates the principal- value 
integral. As usual the real and imaginary parts contribute to the Hall and longitudinal 
conductivities respectively. However, the previous expression would present a singular 
behavior that is an artifact of the n — > limit. This problem is solved by including the 
disorder broadening effects. A formal procedure to obtain the Green function requires a 
self-consistent calculation using the Dyson equation for the self-energy with the magnetic 
and microwave fields, impurity, phonon, and other scattering effects included. A detailed 
calculation of ImG^E) incorporating all these elements is beyond the scope of the 
present work. Instead we choose a gaussian-type expression for the the density of 
states. This expression can be justified within a self- consistent Born calculation that 
incorporates the magnetic field and disorder effects jUHl [3U ES] , hence the density of 
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states for the /i-Landau level is represented as 



ImG,{£) = ,/^exp [-(5-^)7(21^)], (33) 



21* 

with a broadening width given by 

K = ^ (34) 

the parameter (3^ takes into account the difference of the transport scattering time Tj r 
determining the mobility p, from the single-particle lifetime t s . In the case of short- 
range scatterers r tr = r s and (3^ = 1. An expression for (3^, suitable for numerical 
evaluation, that applies for a long-range screened potential is given in reference [T7] ; (3^ 
decreases for higher Landau levels; e.g. /3 ~ 50, j3 50 ~ 10.5. 

The static limit with respect to the external field is obtained taking Q — > in 
Eqs. (13 ip and (J32j) . In what follows results are presented for the microwave-modulation 
(MM) induced longitudinal conductivity, the dark conductivities as well as the MM 
Hall conductivity are quoted in the appendix. Hence the MM induced longitudinal 
conductance is worked out as 

a (MM) = ^-fdSj2J2J2 IrnG„(£)B^ (£,£,) \q%) J x (|A mn |) V mn D, u (q mn )^ (35) 

where the following function has been defined 
d 
d£ 

Notice that a xx contains a contribution o~ xx that depends only on the modulation 
potential, it can be extracted from Eq. (|33|) if the microwave field is switched-off. 

As usual the resistivities are obtained from the expression p xx = a xx J (a1 x + cl y ) 
and p xy = a xy j [a^. x + a^ y ) . The relation a xy 3> a xx holds in general, hence it follows 
that p xx oc a xx , and the longitudinal resistivity follows the same oscillation pattern as 
that of <r rT . 



B®(£,8 V ) = ~{ [f(£ + lu + u E + So) - /(£)] ImG v (£+lu+u E +£o)} .(36) 



4. Results. 

The expression in Eq. (I35J1 can be numerically evaluated after the Fourier components 
V mn of the periodic potential are specified. We shall consider a square lattice potential 
of the form 



V(r) = V 



2irx\ f^wy 

cos I + cos 



(37) 



a 

In our calculations it is assumed that a superlattice is cleaved at the interface 
of an ultraclean GaAs / Al x Ga\- x As heterostructure with high electron mobility, p ~ 
0.5 — 2.5 x 10 7 cm 2 /Vs; the periodic potential has the form given in (|37j) with parameters 
a ~ 20 — 200 nm and Vo = 0.05 meV. The other parameters of the sample are estimated 
as effective electron mass m* = 0.067 m e , fermi energy e^ = lOmeV, electron density 
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n = 3 x 10 n cm~ 2 , and temperature T = IK. For the applied external fields we 
consider: magnetic fields in the range 0.05 — 0.4 Tesla, and microwave radiation with 
frequencies / ~ 10 — 200 Ghz; and field intensity \E^\ ~ 1 — 100 V/cm, that corresponds 
to a microwave power characterized by the dimensionless quantity a ~ ct^E^ 2 / {m* uo z ) 
that varies in the range a ~ 0.01 — 2. The relaxation time r tr in Eq. (|3*3|) is related 
with the zero field electron mobility through \x = eT tr /m*, and (3^ ~ 10.5, a value 
that is justified for large filling factors p, ~ 50 ^7j. A detailed account of the electron 
dynamics requires to distinguish between various time life's; following reference j37j . 
T rad in Eq. (fTUj) is related to the radiative decay width that is interpreted as coherent 
dipole re-radiation of electromagnetic waves by the oscillating 2D electrons excited by 
microwaves, it is given by T rac i = 2ir 2 hne 2 / (3eocm*). In all the examples, except in Fig. 
IH1 we consider the linear regime, the dc-electric field is included only through the Kubo 
formula, hence uje = in Eqs. (|33J) . (jHEJ)- In the case of Fig. IHlthe non-linear dc-electric 
field effects are included using the solution to the classical equations of motion with 
both ac- and dc-electric fields (fTUJ) . 

Plots of the longitudinal and Hall resistivities as a function of the magnetic field 
intensity are displayed in Fig. [T] The total longitudinal resistance shows a strong 
oscillatory behavior with distinctive NRS, this behavior is contrasted with the dark 
contributions that presents only the expected Shubnikov-de-Hass oscillations. The total 
Hall resistance presents a monotonous behavior, yet perceptible microwave induced 
oscillations in the Hall effect can be observed if one considers Ap xy = p xy — p d x a y Tk : , 
see the inset. 

Fig. EJdisplays plots of the total longitudinal conductivity as a function of e = u/u c . 
o xx shows a strong oscillatory behavior, with distinctive negative conductance states. 
The periodicity as well as the number of NCS depend on the intensity of the microwave 
radiation. For weak microwave intensity (a = 0.01), a xx is positive with a moderate 
oscillatory behavior. As the microwave intensity increases (a = 0.1), strong oscillations 
in a xx are observed with minima centered at e ~ 1.2 , e ~ 2.2, e ~ 3.2, and e ~ 4.2; only 
the two last minima correspond to NCS. A further increase in the microwave intensity 
(a = 0.4) yields several NCS. In the region e < 3 the oscillation period is reduced 
to e = |. In general it is observed that a xx vanishes at u>/u> c = j for j integer. The 
oscillations follow a pattern with minima centered at u/u c = j + |(7 — 1) + <5, and maxima 

centered at uj/u c = j + |(/ — 1) —S, where j = 1, 2, 3 , 5 ~ 1/5, and / is the dominant 

multipole that contributes to the conductivity in Eq. ()35j) . For moderate microwave 
power the I = 1 "one photon" stimulated processes dominate, corresponding to what is 
observed for a = 0.1. For a = 0.4 the results can be interpreted as the results of "one 
and two photon" processes (/ = 1 and 1 = 2). To understand the origin of NCS it is 
noticed from Eq. (|33j) that, for small microwave power, a xx is dominated by the / = 
Bessel term, that is always positive. Negative conductance states arise when the / = 1 
and / = terms become comparable: | Jq (|A|) \ 2 B^ ~ | Ji (|A|) \ 2 B^ l \ A simple analysis 
show that this condition is fulfilled for |A| ~ 0.1. Using the expression in Eq. (|2*2*|) . the 
condition to produce NCS can be estimated as IEJ > E th where E th ps 0.1 a T rad /\/8elB- 
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For the parameter used in Fig. |21 E t h ~ 10 V/cm or a t h ~ 0.15, in good agreement with 
the results displayed by the plots. 

Next we consider the dependence of p xx on the lattice parameter a. Plots of p xx 
versus e = uj/uj c for various selections of a are presented in (Fig. EJ). NRS appear only 
in a narrow window of values of a around a*, for which the oscillations amplitude of 
p xx attains its maximum. In the present case: a* ~ 25 nm. This behavior follows from 
the structure of the MM induced longitudinal conductance c^ x M given in (J33j). Taking 
into account the form of Dn V (q mn ) in Eq. (|2*Tj) . the leading dependence of c^ x M on the 
lattice parameter is given approximately by o~^f x M ~ \ exp (— \q\ 2 /2). According to Eq. 
(120 j) q oc I /a, consequently the MM contribution is significant only when a is near to a*, 
that is determined as: a* ~ 7t/b/a/2. This estimation is in very good agreement with 
the numerical results presented in the plots. 

Negative magnetoresistance requires ultra-clean samples, the phenomenon appears 
when the electron mobility exceeds a threshold p t h- Fig. |U displays p xx v.s. uj /uj c 
plots for three selected values of p. For p m 0.5 x I0 7 cm 2 /Vs an almost linear 
behavior p xx is clearly depicted (except in the Shubnikov-deHass region). As the 
electron mobility increases to p « 1.5 x 10 7 cm 2 /Vs, the resistance oscillations are 
clearly observed; however, several NRS only appear when the mobility is increased 
to p ~ 2.5 x 10 7 cm 2 /Vs. Eqs. ()35I36|) contain the main ingredients that explain 
the huge increase observed in the longitudinal conductance (and resistance), when the 
periodically modulated system is irradiated by microwaves and its critical dependence 
on the electron mobility. In the standard expression for the Kubo formula there are 
no Floquet replica contribution, hence u can be set to zero in (|36|). if that is the case 
B"> becomes proportional to the energy derivative of the Fermi distribution, that in the 
T — > limit becomes of the form 8(8 — £p), and the conductivity is positive definite 
depending only on those states lying at the Fermi level. On the other hand, as a result 
of the periodic structure induced by the microwave radiation, B® contains a second 
contribution proportional to the derivative of the density of states: -^ImG v (£ + luf). 
Due to the oscillatory structure of the density of states, this extra contribution takes both 
positive and negative values. According to Eq. (J3l?j) this second term (as compared to 
the first one) is proportional to the electron mobility, hence for sufficiently high mobility 
the new contribution dominates leading to negative conductance states. 

The model is used to test chirality effects induced by the magnetic field. Fig. 
shows a xx v.s. u)/uj c plots for various E w field polarization's with respect to the 
current. It is observed that the amplitude of the a xx oscillations are bigger for linear 
transverse polarization as compared to the longitudinal polarization case. Similarly, the 
oscillation amplitudes are enhanced for negative circular polarization as compared to the 
positive circular polarization results; in particular NCS are observed only for negative 
circular polarization. These results are understood recalling that for negative circular 
polarization and u ~ uj c the electric field rotates in phase with respect to the electron 
cyclotron rotation. 
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Fig. |U1 illustrates the fact that the strong oscillations in o~ xx originates from the 
combined microwave- modulation effects. The dark contribution, Eq. (|38|) . shows 
the expected a xx oc B linear behavior. The contribution arising solely from the 
periodic modulation a^ x can be isolated from Eq. (|3"5)l by switching-off the microwave 
field, the dashed line shows a smooth behavior. The combined microwave-modulation 
contribution is obtained from a^f x M — oi%.\ the dotted line clearly shows that this 
contribution includes the main oscillatory behavior of the full o~ xx , furthermore it is 
the only contribution that becomes negative. Hence, it is concluded that both the 
periodic modulation and the microwave radiation are essential in order to observe the 
strong a xx oscillations and the NCS. 

Next we explore the behavior of the longitudinal conductivity as a function of the 
microwave radiation intensity. As the intensity of the electric microwave field (E w ) 
is increased higher multipole (/) terms needs to be evaluated; in the explored regime 
convergent results are obtained including terms up to the I = 5 multipole. In Fig. [7| 
results are presented for a xx vs E u ; the selected values of e = to Juj c correspond to minima 
or maxima of o~ xx in figure El In general it is observed that for values corresponding 
to maxima i.e. e = j — 5; j = 1,2,3..., a xx remains positive for all microwave field 
intensities. On the other hand, the minima corresponding to e = j + S ; j = 1, 2, 3... 
are related to NCS around a region 1^1 e [10,25] V/cm. Increasing the microwave 
intensity leads to the disappearance of the NCS, except for the first minima (e ~ 1.1). 

The non-linear regime with respect to the dc-external field can also be explored 
within the present formalism. The effect is included using the solution to the classical 
equations of motion with both ac- and dc-electric fields, Eq. (JIDj) . A possible connection 
between the observed ZRS in GaAs / Al x Ga\- x As heterostructures [U El El Ej and the 
the predicted NRSpm HH H21 US dll HH1 CHI [HI was put forward by Andreev etal. 
[T5] . noting that a general analysis of Maxwell equations shows that NRS induces an 
instability that drives the system into a ZRS. This mechanism requires the longitudinal 
current j xx as a function of E^ c to have a single minimum, the system instability will 
evolve to the value Ed c in which j xx cancel. Returning to the irradiated superlattice case, 
in Fig. [HI it is observed that in general the j xx vs. E^ plot has an oscillatory behavior, 
with more than one minima. Hence the conditions of the Andreeev mechanism do 
not apply. Consequently, negative conductance states may be probably observed in 
2-dimensional superlattices, when exposed to both magnetic and microwave fields. 

5. Conclusions. 

We have considered a model to describe the microwave photoconductivity of a 2DES 
in the presence of a magnetic field, and a 2D periodic modulation. We presented a 
thoroughly discussion of the method that takes into account the Landau and microwave 
contributions in a non-perturbative exact way, the periodic potential effects are treated 
perturbatively. The formalism exploits the symmetries of the problem: the exact 
solution of the Landau-microwave dynamics ((JJ) is obtained in terms of the electric- 
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magnetic generators (J5J) as well as the solutions to the classical equations of motion 
(jHJ). The spectrum and Floquet modes are explicitly worked out (fTTj) . In our model, 
the Landau-Floquet states act coherently with respect to the oscillating field of the 
superlattice potential, that in turn induces transitions between these levels. Based on 
this formalism, a Kubo-like formula is provided (|33|) that consistently takes into account 
the oscillatory Floquet structure of the problem. 

It is found that both a xx and p xx exhibit strong oscillations governed by e = uj/uj c . 
The oscillations follow a pattern with minima centered at uj/oj c — j + |(Z — 1) + S, and 

maxima centered at uo/uj c = j + \{l — 1) — S, where j = 1,2,3 , 5 ~ 1/5 and I is 

the dominant multipole contribution. NRS develop for sufficiently strong microwave 
power (Fig. EJ), in a narrow window of values of the lattice parameter (a ~ Ib) (Fig. |3J), 
and for high electron mobility samples (Fig. HJ). The explanation for the NRS can be 
traced down to Eqs. (j33|) and (jHUj) ; the longitudinal photoconductivity contains a new 
contribution proportional to the derivative of the density of states: -^ImG v {£ + lu). 
Due to the oscillatory structure of the density of states this extra contribution takes 
both positive and negative values. This term is proportional to the electron mobility, 
hence for sufficiently high mobility the new contribution dominates leading to negative 
conductivity states. Unlike the semiclassical origin of magnetoresistance oscillations 
observed in an antidot array for commensurate values of the ratio R c /a J2H], these 
conductance oscillations have a quantum origin and would only appear in a narrow 
window of values of a, around a ~ Ib- 

In conclusion, it is proposed that the combined effect of: periodic modulation, 
perpendicular magnetic field, plus microwave irradiation of 2DES give rise to interesting 
oscillatory conductance phenomena, with the possible development of NCS and NRS. 
One should stress that according to our results, the production of NRS requires ultra- 
clean samples with electron mobilities of order fi « 2.5 x 10 7 cm 2 /Vs (see Fig. HJ). 
The electron mobilities in the fabricated arrays of periodic scatterers so far |2"H l2"K] are 
p fh 2.5 x 10 6 cm 2 /Vs, consequently an increase on the electron mobilities of these kind 
of experimental setups by an order of magnitude would be required in order to observe 
the phenomena described in this work. 

5.1. Appendix A: Dark and Hall conductivities. 

In section (jSJ) it was explained in detail the method to obtain the final expression for 
the modulation-microwave induced conductance Eq. (|35|) . Working along a similar 
procedure the expression for the remaining conductivities are worked from equations 
(j3*Tj) and (JH2J)- First we quote the longitudinal dark conductance 

,2, ,2 



D e _^c 



A' 

whereas the Hall conductance reads 
n e 2 u 2 c 



^2fi fdSImG^S) ^ImG^S + Uc) , (3* 



"■:.. = —r E A* / d£ Im G, {£) [f {£, -u c )-f {£)] P - l , (39) 
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where P indicates the principal- value integral. The final result for the MM assisted 
longitudinal conductivity was quoted in Eq. (|35|). The MM induced Hal conductivity 
is calculated to give 



7rh , 

fiu I ran 

were the function T l mn is defined as 

2 



2 



a (MM) = ^| I d£ J2 J2 J2 ImG » (£ ) [/ (£0 - / (£)} T L Jl (lA^I) V mn D, u (q mn ) ", (40) 



(g^) 2 + tf } 



T in = ^c " " 5 • ( 41 ) 

c (f + o;Z-^) |(£ + cu/-^) 2 -cu 2 | 2 V ; 

5.#. Appendix B: Microwave- driven distribution function. 

Within the time relaxation approximation the Boltzmann equation can be written as 

dt dp Ttr Ti n 

where fp is the Fermi-Dirac distribution and we distinguish between the elastic rate 
r t 7 and inelastic or energy relaxation rate t~^ . As it already mentioned, we assume the 
validity of the following conditions: r^ <C r tr <^ r in , and certainly the inelastic processes 
can be safely ignored. Furthermore, due to the ac-electric field (jHJ), the L.H. S. of the 
previous equation is estimated to be of order //r w ; hence, in a first approximation the 
elastic scattering contribution can also be neglected. The resulting Vlasov equation 
has the exact solution /(p, t) = fp (p — m*v(t)), where the velocity v(£) = (»)i, Ci ) 
solves exactly the same classical equations of motion that follow from (JHJ), and the initial 
condition is selected as, /->/j?as the external electric field is switched-off. In particular 
it is verified that m*|v(t)| 2 /2 = £ rad coincides with the Floquet energy shift produced by 
the microwave radiation (JTTJ). The steady-state distribution, evaluated at the Landau 
energy £ = £ji , is obtained by averaging fp (p — m*v(t)) over the oscillatory period 



(fp) = y ]" fF fc 0) + £ rad + 2 coscu c t y/e^£ rad \ dt. 



(43) 



In general it is verified that £ ra d <^. £ji , thus expanding to first order one finds 
(f F ) w f F ( £\i + £ ra d) = fp(£fj,)- Hence, it is verified that a rapid relaxation of the 
Fermi distribution to the quasi-energy states is a reasonable assumption. The arguments 
presented in this appendix have been introduced by Mikhailov 37^ in order to explore 
the possibility that the microwave radiation leads to a population inversion; however, it 
is concluded that it would require a rather high microwave intensity £ ra d > £p. 
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Figure 1. Magnetoresistivity p xx versus B without radiation and modulation (dotted 
line) and under microwave radiation plus periodic modulation (solid line). The dashed 
line corresponds to the total Hall resistance p xy , whereas the inset shows the results for 
the subtracted Hall resistance Ap xy — p xy — p x °? k . The microwave polarization is linear 
transverse (with respect to the current), with / = 100 Ghz and a = ce^E^] 2 / (m* ui 3 ) 
— 0.4. The other parameters are selected as follows: a — 25nm, Vo — 0.05 meV, 
m* = 0.067 m e , n « 2.5 x 10 7 cm 2 /Vs, e F = 10 meV, T=1K. 
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Figure 2. Longitudinal conductivity versus e = lo/uj c for three values of the microwave 
power intensity: a = 0.01 dashed line, a = 0.1 dotted line, and a = 0.4 continuos line. 
The other parameters have the same values as in Fig. ^ 
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Figure 3. Longitudinal resistivity versus e = w/w c for three values of the lattice 
parameter: a = 22 nm dotted line, a = 25 nm continuos line, and a = 35 nm dashed 
line. The microwave intensity is a — 0.4 and the other parameters have the same 
values as in Fig. Q] 
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Figure 4. Longitudinal resistivity as a function e = uj/lu c for three values of the 
electron mobility: [i — 0.5 x 10 7 cm 2 /Vs dotted line, /U = 1.5 X 10 7 err? jV s dashed 
line, and /i = 2.5 x 10 7 cm 2 /Vs continuos line. The microwave power is given by 
a = 0.4, the other parameters have the same values as in Fig. ^ 
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Figure 5. Longitudinal conductance vs. e — ui/ll> c for various microwave E^— field 
polarization's with respect to the current. In figure (a) the continuos and dotted lines 
correspond to linear transverse and longitudinal polarization's respectively. Figure 
(b) shows results for circular polarization's: left-hand (continuos line) and right-hand 
(dotted line) . a = 4 and the values of the other parameters are the same as in figure 

Fig.m 
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Figure 6. Contributions to the total a xx (continuos line) vs. e = lo/lu c : dark 
contribution a xx (dashed-dotted line) from Eq. i|38|l , periodic modulation contribution 
a^f x (dashed line) obtained from Eq. Q35J I when the microwave field is switched-off, 
microwave-modulation contribution obtained as <J XX — cr xx 
parameters are the same as in figure Fig. ^ 
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Figure 7. Longitudinal conductivity as a function of the microwave ac-electric field 
for various values of e = ui/uj c : e = 1.1 continuos line, e = 1.9 dotted line, e = 2.15 
dashed line, and e = 4.15 dashed-dotted line. The other parameters have the same 
values as in Fig. Q] 
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Figure 8. Current-voltage characteristics for the irradiated sample for various values 
of e = uj/lo c : e = 1.1 continuos line, e = 2.1 dashed line, and e = 3.15 dotted line. The 
microwave power is given by a = 0.4, the other parameters have the same values as in 

Fig.m 



